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where u = (u 1 , M m )* is a vector valued function and a a P(x) are (n X n) matri- 
ces whose elements a"P(x) are at least uniformly bounded measurable real-valued 
functions such that 



for any combination of a,/3,i, and j. We assume we have two non-empty, open, dis- 
joint, and bounded sets, Q and Q, in R n , and add a set T e of small measure to form the 
domain f2 e . Then we show that as e — ¥ + , the Dirichlet eigenvalues corresponding to 
the family of domains {f2 e } e >o converge to the Dirichlet eigenvalues corresponding 
to tto = Q U Q. Moreover, our rate of convergence is independent of the eigenvalues. 
In this paper, we consider the Lame system, systems which satisfy a strong ellipticity 
condition, and systems which satisfy a Legendre-Hadamard ellipticity condition. 
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1 Introduction 

There is a great deal of work studying eigenvalues for elliptic equations, but 
there seems to be less work on eigenvalues for elliptic systems. Much of the work 
on equations requires estimates for solutions that do not hold for systems. In this 
paper, we consider the behavior of eigenvalues for elliptic systems in singularly 
perturbed domains. We give a simple characterization of the families of domains 
that we can study and this class includes families such as dumbbell domains formed 
by connecting two domains by a thin tube. We show that as the measure of the 
perturbation shrinks away, the convergence of the eigenvalues is obtained. We also 
provide a rate of convergence, which is independent of any eigenvalue. We make 
no assumption on the smoothness of the coefficients and only mild assumptions 
on the boundary of the domain. 
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Studying solutions of elliptic boundary value problems with Dirichlet or Neu- 
mann boundary conditions on domains which can be approximated by solutions 
on simpler domains has been an interest for many years, and is still ongoing. The 
motivation to study such problems is that it is easier to study the spectra on 
sets with a reduced dimensionality. One may approximate the spectra on these 
"fattened" sets with the spectra on the "thinner" sets. Some applications include 
studying quantum wires, free-electron theory of conjugated molecules, and pho- 
tonic crystals. For a complete description, see the work of Kuchment [55]. Recent 
work by Exner and Post |15| study the Neumann Laplacian on manifolds with 
thin tubes which is related to the theory of quantum graphs. The Fireman's Pole 
problem consists of approximating the resolvents of a bounded set in M 3 by the re- 
solvents of this set with a cylinder removed. For a complete description, see Rauch 
and Taylor [29]. A classic paper by Babuska and Vyborny [5] shows continuity of 
Dirichlet eigenvalues for elliptic equations under a regular variation of the domain, 
but gives no rates of convergence. Dancer |12) considers how perturbing the 
domain affects the number of positive solutions for nonlinear equations with Dirich- 
let boundary conditions and includes the case where solutions are eigenfunctions 
for the Laplacian. Davies [14] and Pang [28] study the approximation of Dirichlet 
eigenvalues and corresponding eigenfunctions in a domain f2 by eigenvalues and 
eigenfunctions in sets of the form R(e) = {x 6 ft : dist(a;, dft) > e}. They each 
give rates of convergence and their estimates include the case when the domain is 
irregularly shaped. The work of Brown, Hislop, and Martinez [7] provides upper 
and lower bounds on the splitting between the first two Dirichlet eigenvalues in a 
symmetric dumbbell region with a straight tube. Chavel and Feldman [9] examine 
eigenvalues on a compact manifold with a small handle and Dirichlet conditions on 
the ends of the handle. The work of Anne and Colbois [1] examines the behavior 
of eigenvalues of the Laplacian on p-forms under a singular perturbation obtained 
by adding a thin handle to a compact manifold, but requires more regularity on 
the eigenfunctions than holds in our setting. 

More recent work for Dirichlet conditions includes work by Daners [T3] , which 
shows convergence of solutions to elliptic equations on sequences of domains. These 
domains 0„ converge to a limit domain Q in the sense of sequences u n 6 i/g(f2 n ) 
converging to a function u S Hq(£1). Also, Burenkov and Lamberti [5] prove sharp 
spectral stability estimates for higher-order elliptic operators on domains in certain 
Holder classes in terms of the Lebesgue measure of the symmetric difference of the 
different domains. Kozlov [25] obtains asymptotics of Dirichlet eigenvalues for 
domains in W for n > 2 using Hadamard's formula. Grieser and Jerison [20] also 
give asymptotics for Dirichlet eigenvalues and eigenfunctions, but only on plane 
domains. 

We note here that the results for Neumann eigenvalues may be different than 
those for Dirichlet eigenvalues. In fact, a classic example of Courant and Hilbert 
[10] shows that the Neumann eigenvalues may not vary continuously as the domain 
varies. Their example is constructed by taking the unit square in R 2 and attaching 
a thin handle with a proportional square attached to the other end. They show 
that if {A^} and {A°} are the Neumann eigenvalues of —A in increasing order 
including multiplicities with respect to the unit square and the perturbed square, 
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then A| — > as e — > 0, but A° > 0. This example shows that one needs additional 
regularity in order to achieve convergence. Furthermore, Arrieta, Hale, and Han [3] 
show that for this type of domain, A^ — > A^ 1 _ 1 , as e — » for m > 3. Another work 
of Arrieta [1] gives rates of convergence for eigenvalues of the Neumann Laplacian 
on a dumbbell domain in M 2 when the tube is more general. Jimbo and Morita 
|23] study the first N eigenvalues of the Neumann Laplacian in N disjoint domains 
connected by thin tubes. They show that the first N eigenvalues approach zero 
and the (JV + l)st eigenvalue is uniformly bounded away from zero. If D\ and D2 
are two disjoint domains, then for {c^} = {/i/}U{Aj}, where {/i/} are the Neumann 

eigenvalues of — A in D = D\ UZ?2 and {Xj} are the Dirichlet eigenvalues of -^pr in 
(—1,1), Jimbo [22] gives a rate of convergence on the difference a\. — erf. This work 
was generalized to more classes of domains in a more recent work by Jimbo and 
Kosugi [21]. Also, Brown, Hislop, and Martinez [3] show that if <ik £ {^i}\{Xj} 
then 

-1 

2 

n = 2 

|cr fe - erg I < Ce^ n>3. 

Here, we aim to provide an outline of the proof. In section 2, we give several 
definitions and describe the family of domains for which we can prove the conver- 
gence of eigenvalues. We also describe the well-known construction of eigenvalues 
and state our main result. In section 3, we give Theorem 13 . 1 1 from Giaquinta and 
Modica [T7], [TB] which uses a technique introduced by Gehring [TB]. We also 
prove a Caccioppoli type estimate for eigenfunctions in Theorem 13.41 and use this 
along with Theorem 13.11 to obtain a reverse Holder inequality given in Theorem 
13.51 This gives i p -integrability for the gradient of the eigenfunctions for p > 2. In 
section 4, we are able to bound these LP norms by a constant in Proposition 14.21 
The proof uses the reverse Holder inequality as the key ingredient. This estimate 
is then used to prove Lemma |4~21 and Proposition 23] which are used to satisfy the 
first part of a well-known theorem from Anne [2] given in Lemma |4~T1 The second 
part of Lemma |4 . 1 1 follows from the first part along with the above estimates, thus 
giving Corollary 14. 11 The main result follows from this corollary. As a by-product 
of our research, we give a simple proof of Shi and Wright's (30] L p -estimates for 
the gradient of the Lame system as well as other elliptic systems. Many of the 
results first appeared in the author's Ph.D. dissertation |31| . 
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2 Preliminaries and Main Result 

We give conditions on a family of domains Q e that allow us to prove the 
convergence of eigenvalues. We let f2 and 11 in R™ be two non-empty, open, 
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disjoint, and bounded sets. We let e\ > (which will be chosen small later), 
and then let {T £ }q <£ < £i be a family of open sets such that 

X~ c T £ if e < e 

and if \T e \ denotes the Lebesgue measure of T £ , then 

\T e \ < Ce d (2.1) 

where C and < d < n are independent of e. Fix two points p\ and pi on <957 and 
<957, respectively. For each e, let B £ and B £ be two balls of radius e in 1" centered 
at pi and p2, respectively. The connections from T e to 57 and 51 will be contained 
in B £ and B £ , so that T £ n 57 = and T £ nUc B| where .Bj is the concentric ball 
to B £ of radius |. Also, suppose a similar condition for 57 and B £ . Then for any e, 
define 57 e to be the set 57 U 51 U T £ , which we assume to be open and connected, and 
57o = 57 U 57. So, if our family is the family of dumbbell domains, you may think 
of T £ as a "tube" connecting each of the two domains. We now have the family of 
domains {57 e }o< £ < £l . 

Next, we give a condition on the boundary of 57 e . If B r is any ball of radius r 
satisfying B r n £l £ ^ 0, then 

\B 2r n 57^| > C* r™ (2.2) 

where Co is a constant independent of r and e. This eliminates domains with 
"cracks" and "in-cusps," and will be used to help show the Caccioppoli inequality 
in Theorem 13.41 for the case when we are close to the boundary. 

Throughout this paper we use the convention of summing over repeated indices, 
where i and j will run from 1 to n and a, /3, and 7 will run from 1 to m. We 
let a?j (x) be bounded, measurable, real- valued functions on l n which satisfy the 
symmetry condition 

a ij( x ) = a ji{ x )i h3 = l,2,...,n, a,P = 1,2, ...,m. 

We let L 2 (n £ ) denote the space of square integrable functions taking values in 
WL m and 7Jq(57 £ ) denotes the Sobolev space of vector- valued functions having one 
derivative in L 2 (Q £ ) and which vanish on the boundary. We use u" to denote the 
partial derivative 

Let r] £ G C£°(K") be a cutoff function so that r] e = in T £ , 7? e = 1 in 57 \(B e U 
£> e ), |V»7e| < g j an d < ?7 £ < 1, where C n only depends on n. We emphasize 
that B £ , B £ , and ?7 6 depend on the parameter e. With these assumptions and 
definitions, we have that for any u G i?Q(i7 e ), will be in ^(f^o). 

We now introduce the notion of an eigenvalue and corresponding eigenvector. 
We say that the number a is a Dirichlet eigenvalue of L with Dirichlet eigenfunction 
u G H$(n), if u ^ and 

/ aff(x)uf {x)<j)?(x) dx = a \ v? '{x)<jp ' (x) dx, for any € Hq(SI). 

Jn Jn 

(2.3) 
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We say that L satisfies the Legendre-Hadamard condition if there exists 9 > so 
that 

a?f(x)USpM, >0|£| 2 H 2 , ^« n , a.e. x G £ . (2.4) 

If we define the norm on matrices A = Ay G K mxn a s |A| 2 = | A i3 ■ | 2 , and L 

satisfies the Legendre-Hadamard condition with continuous coefficients in f2, then 
it is well-known that for any u G Hq(O), we have Garding's inequality [521 p. 347] 

Ci / \Vu\ 2 dx< [ a™f(x)u?(x)u?(x)dx + C 2 [ \u\ 2 dx. (2.5) 
Jn Jn Jn 

L is said to satisfy a strong Legendre condition or a strong ellipticity condition if 
there exists > so that 

«<? > ^l 2 , eeK mxn , (2.6) 

We introduce the Lame system as Lu = — div£(u), where £(tt) denotes the 
stress tensor defined by 

C/(u):=a^uf (2.7) 
which is defined in terms of the Lame moduli and n(x) by 

°£f 0*0 = v{x)5 ia 8jp + n{x)8ij5 a p + fi(x)SipS ja7 (2.8) 

where v{x) and /i(a;) are both assumed to be bounded and measurable. Also, 
define the strain tensor k(u) as 

Note that for the Lame system, m = n and the Lame parameters v(x) and n(x) 
given in (|2.8|) satisfy the conditions 

v(x) > fi{x) > 5 > 0. (2.10) 

With these assumptions, the Lame system satisfies the ellipticity condition 

> t \k(u)\ 2 , u G H^(Q £ ) (2.11) 

where r = 2(5. With Korn's 1st Inequality, it is easy to see that for the Lame 
system, we have 

T - [ \Vu\ 2 dy<[ affuf^dy, u G Hq(Q e ). 
1 Jn c Jn E 

Thus, if u satisfies either the ellipticity condition (I2.6p . (|2.1ip . or (|2.4D with con- 
tinuous coefficients in f2, then we have Garding's inequality (|2.5| . 
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The well-known construction of eigenvalues and eigenfunctions for scalar func- 
tions (which is the same for vector-valued functions) is taken from Gilbarg and 
Trudinger [19, p. 212]. If we define the bilinear form on Hq(£1 £ ) x H^(fl e ) as 

B E (u,v):= [ affufv? dx (2.12) 
Jn e 

and define the Rayleigh quotient R £ as 

B E (u,u) 

R z\ u ) ■= Tj — jTa ( 213 ) 

for u ^ 0, then we can construct an increasing sequence of eigenvalues, listed 
according to multiplicity, {(Jk}^ 1 such that for each corresponding eigenfunction 
Ufc £ Hq(£1 £ ), we have 

min R s (w) = R E (u k ) = o k (2.14) 
ioe{iii,...,Ufc_i}- L 

and 

ll«*IU»(n.) = 1 (2-15) 

for any k. Furthermore, each eigenspace is finite-dimensional and the constructed 
set of eigenfunctions forms an orthonormal basis in L 2 (Q, e ). 
We now state the main result. 

Theorem 2.1. Let 

satisfy one of the following: 

1. L has uniformly bounded coefficients and satisfies either the ellipticity con- 
dition A2. 6]) or the ellipticity condition i2.11]) . 

2. L has continuous coefficients and satisfies the ellipticity condition \2.J$ . 

Also assume {o~k}kLi an d {"felfeLi are Dirichlet eigenvalues of L with respect 
to flo and fl e in increasing order numbered according to multiplicity. Then for 
each J € N, we have the following estimate: 

\a^-a° J \<Ce a 

for < £ < £o(J), where £o(J) depends on the multiplicity of a®. Moreover, the 
rate a > is independent of any eigenvalue and C only depends on the eigenvalue 
a j and the distance from a , to nearby eigenvalues. 
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3 A Reverse Holder Inequality 



if 



T \f(y)\ dy is defined to be the average of / on E, then recall that the 
Je 



maximal function is defined for / G Lj oc (M. n ) to be 

M(f)(x) :=sup-/ \f(y)\dy 

r>0 J B r {x) 

where B r (x) is a ball of radius r centered at x. Also, define Mn{f){x) to be 



M R {f)[x) := sup 4 \f(y)\ dy. 

R>r>0J B r (x) 

We will need the following theorem from Giaquinta |17[ p. 122], which uses the 
technique introduced by Gehring |16| . and refined by Giaquinta and Modica [18J. 

Theorem 3.1. Let r > q > 1, and Qr be a cube in R" with sidelength R centered 
at 0. Also, define d{x) = dist(x,dQn). If f and g are measurable junctions such 
that f G L t (Qr), g G L q {Qii), f = g = outside Qr, and with the added condition 
that 

M^{\g\i){x) < bM«(g)(x) + M(\f\«) + aM(\g\«)(x) 

2 

for almost every x in Qr where b > and < a < 1, then g G L p (Qr), for 
p G [q,q + e) and 



J \g\ p {y) dy) <c 



\g\ 9 (y)dy) +14 \f\ p (y)dy 



(3.1) 



where e and C depend on b, q, n, a and r. 

The conclusion of this theorem is known as a reverse Holder inequality. To 
show that the gradient of eigenfunctions satisfy this inequality, we will need to 
prove a Caccioppoli inequality. However, to show this Caccioppoli inequality, we 
first need the following two well-known inequalities taken from Hebey [2lT p. 44] 
and Oleinik [27l p. 27]: 

Theorem 3.2. Sobolev-Poincare Inequality Let 1 < p < n and ~ = ~ — i. 

Also, let B r be any ball of radius r with u G W 1,p (B r ). Then, for S contained in 
B r with \S\ > c r n , 

\u(x) - u s \ q dx < C (J \\7u\ p (x) dxj (3.2) 

where us = j^y J s u dy, for some constant C(n,p,co), independent of u. 
Theorem 3.3. Korn's Inequality on a Ball If u E H 1 (B r ) then 



|Vu||i a(Br) < C ( Mu)\\l 2{Br) + ±\\u\\l 2{Br) ) (3.3) 



where C only depends on n. 
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We now state and prove a Caccioppoli inequality for eigenf unctions: 



Theorem 3.4. Let u be an eigenfunction with eigenvalue a associated to the 
operator L satisfying either \2. 6}) or i2.11\) with uniformly bounded coefficients or 
associated to \2.1$ with continuous coefficients. Extending u to be outside Q £ , 
there exists r > so that if > r > 0, x £ W 1 , we have 



£ |Vm| 2 dy <d U |Vtt|^ dy^j 

+ C 2 \a\-f \u\ 2 dy + C 3 -f |Vu| 2 dy (3.4) 

J B 2r J B 2r 

where B r is a ball with radius r centered at x, C3 < 1, and Ci > only depends 
on M — maxi JiQj ^ ||oiy ||i,=«<(o e ), n, m, 9, t, and Co. Furthermore, if L satisfies 
either \2. 61) or i2.11)) with uniformly bounded coefficients, then the inequality holds 
for any r > 0. 

Proof. First, choose a ball B r and define a cutoff function v £ C£°(R") to be so 
that v = 1 in B r , v = outside B<z r , \Wv\ < and < v < 1, where C„ only 
depends on n. Below, we will find an appropriate constant vector p £ R m , so that 
v 2 (u — p) £ Hq(£1 £ ). By the weak formulation (|2.3[) . we have 

n e ./n e 
Then, performing the differentiations, we get 

a°f <[2z^(u - + i/V] dy = a [ v?v 2 {u - P y dy. (3.5) 



From this point, the argument depends on the ellipticity condition. We have 3 
cases. 

case 1: L satisfies the strong ellipticity condition 112. 6\) . 

Using (|2.6p and properties of v, we obtain the inequality 

v 2 afj 1 u" Uj dy < [ 2M— —v\ Vtt| \u — p\ dy + ( \a\\u\\u — p\dy 

B 2r B 2r ' J B 2r 

which, for any constant u) > 0, then leads to 

+ C\a\ f \u\ 2 dy (3.6) 

JB 2r 
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where C depends on M and C n . Then choosing uj = 9 in Q3.6P gives 

°-( v 2 \Vu\ 2 dy < f \u-p\ 2 dy + C\a\ [ \u\ 2 dy. 

L JB 2r 0r J B 2r JB 2r 

Then, multiplying both sides by | and using that v = 1 on B r gives 

\Vu\ 2 dy<^[ \ u -p\ 2 dy+ 2 -^ [ \u\ 2 dy. (3.7) 



Now. 



for the term / \u — p\ 2 dy, we must consider two subcases. 



ibcase A 



If Bir C f2 e , then let p a — 4 u a dy. Our condition on the support of v 

Jb 2t 

implies v 2 {u — p) G Hq(VL e ). So, setting q = 2 and S = B^r in the Sobolev- 
Poincare Inequality (|3.2p . we obtain 

/ \u-p\ 2 dy<c( f \Vu\^ dy) " . 
Using this estimate with Q3.7P gives 

re + 2 

j B \^\ 2 dy<^J g \Vu\^ dy) " +CW\J b M 2 dy. 
Now, dividing through by r n gives the desired result with C3 = 0. 

subcase B 



If i?2r H 7^ 0, then set p = 0, which, again, guarantees that ^ 2 (m — p) £ 
-£Tq (f2 e ). So setting g = 2 and S 1 = i?4 r n f2 e in the Sobolev-Poincare Inequality 
we have by our assumption on VL c e (|2.2p that 



f \u-p\ 2 dy<c(f \\7u\^2 dy 
From (|3.7p . we obtain 

Tl+2 

/ |Vu| 2 dy<^( f |Vu|3* dy ) " + / M 2 dy. 

JS r r \J B ir / •'■B.ir 

A simple covering argument gives the estimate with B^ r replaced with Bi r . 
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2: L satisfies the ellipticity condition i2.11\) . 



From (|2.11[) and (|3 . 6f) . we have 
r\ K (u)\ 2 dy< / l - — - 

B r JB 2r z 



d V+^l \u - p\ 2 dy + C\a\ [ \u\ 2 dy. 
ur 2 J B2r Jb 2t 



Also, by Korn's inequality (|3.3p , we have 

IjT Nu\ 2 dy-^J^ \u~p\ 2 dy<J^ t\k{u)\ 2 dy. 

This implies 

/ Nu\ 2 dy<^[ \Vu\ 2 dy + c(^ + ±) f \u-p\ 2 dy 

JB r ZT J Bo,, \UJTT T J J g 



C\a 



B 2r 

u\ 2 dy. 



B-2 



This again leads to two subcases. We must choose p appropriately and use the 
Sobolev-Poincare inequality (|3.2|) as in case 1. Then, by taking lu sufficiently small, 
we obtain the desired result. 



case 3: L satisfies the Legendre-Hadamard condition (|2.4[) with continuous coef 
cients in fL. 



We note that it suffices to study when u <G C%°(Sl £ ) and first consider when the 
coefficients arc constant. We rewrite the left side of (13.51) as 



<%? {{u-pYvUiu-pfv^dy 



+ J a"f[ui/juf(u - p) p - Viv[u - p) a u^ - ViVj(u - p) a {u - p) p ] dy. 
This implies 

afdu-pYvUiu-pY^dy 

B 2r 

<C |Vf||V((u - p)v)\\u - p\ + \u - p| 2 |V^| 2 + \cr\\u\\u - p\ dy. 

JB 2r 

We note that we may use the Fourier transform to get a lower bound for the left 
side to achieve the estimate 

/ \V((u-pM 2 dy<^ [ \u~p\ 2 dy + C\a\ [ \u\ 2 dy. 

J B 2t r J B 2r J B 2r 
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This implies the estimate 

f \ Vu \ 2 dy<^[ \u-p\ 2 dy + C\a\ f \u\ 2 dy. (3.8) 

J B r r J B 2r J B 2r 

So, again, if we employ the Sobolev-Poincare inequality (|3.2p . we get the desired 
result in the case of constant coefficients. If the coefficients are continuous and non- 
constant, then we freeze the coefficients at x. That is, from the weak formulation 
(|2.3j) . we have 

/ a?f(x)u?((u p)^ dy + / (aff - a?f (x))u? ((« - p)^ dy 

= a [ u' f ((u- p)v 2 y dy. (3.9) 

So, if we define the modulus of continuity to be 

M(x ,R) = max | a°f (y) - a°f (x )\ 

y£B R (x ) 

then we have that 

j B K?-af{x))u?{(u-p)v 2 f 3 dy 

<M(x,2r) v 2 \Vu\ 2 dy + 2M(x 7 2r) v\ Vv\ | Vit| \u - p\ dy 
<C(M(x,2r)+M(x,2r) 2 ) [ \\7u\ 2 dy+^[ \u-p\ 2 dy. 

JB 2r r JB 2r 

Also, by the uniform continuity of the coefficients on f2 e , for any c < 1, there exists 
ro depending on c, so that if C(xq,R) = C(M(xo,2R) + M{xq 1 2R) 2 ) and r < r , 
then 

C(xo,r) < c 

for all xo 6 Sl £ . So, now moving the second term on the left side of (13 ,9|) to the 
right and using the constant coefficient case (|3.8|) . we obtain that for any c < 1, 
there exists ro so that if r < ro , 

f \Vu\ 2 dy<^[ \u - p\ 2 dy + C\a\ [ \u\ 2 dy + c [ |Vm| 2 dy. 

J B r r J B 2r J B 2v J B 2r 

We again choose p appropriately and apply the Sobolev-Poincare inequality (|3.2p 
to get the desired result. □ 

As stated earlier, our proof of Theorem 12.11 relies on the gradient of an eigen- 
function satisfying the reverse Holder inequality, as in our next theorem. 



11 



Theorem 3.5. There exists e\ > so that if u is an eigenfunction with eigenvalue 
a, then 



/ 



|Vw| p dy < C 



|Vu| 2 dy 



\u\ p dy 



(3.10) 



where 2 < p < 2 + e l7 and ei and C are independent of e and any eigenvalue. 

Proof. Now if u is an eigenfunction with eigenvalue a, we have u £ H^fl^), and 
thus we may employ the Sobolev inequality to get that \u\ £ L r (Vl e ) for some r > 2. 
If L satisfies either (|2.6I) or (|2.1ip with uniformly bounded coefficients, then we 
may choose a cube Qr, centered at 0, with sidelength R such that fl e C Qr, 

uniformly in e. and set g = |Vw|^, / = (C 3 \a\)^\u\^ , q = and u = 
outside f2 e , we may conclude by (|3.4[) and (|3.1[) that 



2rep 

Vtt dy 



< C 



\Vu\ 2 dy 



a h+ 2 



where < p < s±£ ^_ £; w hich, from Theorem 13.41 is independent of e and 
any eigenvalue. So, setting p = we have the result. If L satisfies (|2.4[) with 



continuous coefficients, then since we only have Theorem 13.41 true for small r, we 
must cover f2 e with a fixed number of cubes and apply ()3.1[) to each cube to obtain 
the result. □ 



4 Stability of Eigenvalues 

From this point, let <r| be the kth eigenvalue with respect to fL, and 0| be 
its corresponding eigenfunction with (f) £ k — outside J7 e for e > 0. We also fix an 
eigenvalue o~j with multiplicity mj where cr°j_i < o~j if J > 2. We will consider 
the family {oj} as e > tends to 0. We begin with the following proposition 
taken from Anne d p. 2595-2596]. 

Lemma 4.1. Let (q,T>) be a closed non-negative quadratic form with form domain 
T> in the Hilbert space (~H, (•, •)). Define the associated norm \\f\\i = \\f\\^i +?(/), 
and the spectral projector 11/ for any interval I = (a, fj) for which the boundary 
does not meet the spectrum. 

1. Suppose f £ T> and X E I satisfy 

W, g) -Mf,g)\<s\\f \\\\9h gev. 

Then there exists a constant C > 0, which depends on I, such that if a is 
less than the distance of a or fj to the spectrum of q, 

\\U I (f)-f\\ 1 = \\U I .(f)\\ 1 <—\\f\\. 
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2. Suppose the spectral space E(I) has dimension m and /i,...,/ m is an or- 
thonormal family which satisfies 

\\n.i'Ui)h<& ./ i 

Also let E be the space spanned by the fj 's. Then, 

dist(E(I),E) < CS 

where the distance is measured as the distance between the two orthogonal 
projectors. 

This lemma will give us the results we need for the convergence of eigenvalues. We 
will prove estimates on eigenfunctions using the reverse Holder inequality (|3.10p . 
which will allow us to use this lemma. We start with the following proposition 
which follows immediately from the construction of eigenvalues. 

Proposition 4.1. We have for any e > 0, and any k £ N, 

°i < 4- (4.i) 

This proposition gives us the easy half of the inequality in our theorem. To prove 
the second half of the inequality, we will need a few items. 

Proposition 4.2. For any e > 0, and k > 1, if (f> = (/)%, then we have 

\V(f>f dy<C (4.2) 

where p > 2 is from (|3.10p . andC depends on |fio| andn, with order O ^|cr°| P+ i P ' 

( qp + 2(p-q) \ 
|oj!| 2 i ] for n = 2 where 2 — £ < q < 2 for small £. Fur- 
thermore, p and C are independent of e and ifn — 2, C blows up as q — > 2. 

Proof. Now, from (I3.10|) . we have 



f \V<pfdy<C \n e \^(j |V0| 2 dy) ( f Wdy) 



(4.3) 



where p > 2 is from (|3.10p . Recall that by Garding's inequality (|2.5p and since 4> 
is an eigenfunction, we have 

Cx [ |V0| 2 dy < [ aftfrf dy + C 2 [ |0| 2 dy 
<C(l + \a%\) f H 2 dy 

<C(l+\o%\), (4.4) 
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the last line owing to the normalization of the eigenfunctions. Next, we will con- 
sider n > 3 and estimate 



W dy. 

Using Sobolev's inequality and (|4.4p . we have 



-2 dy 



Also, by Holder's inequality, we have 



< C U \V4>\ 2 dy 
<C{l + \a%\i). 



dy) <[ W dy 

la. 



I"- 2 dy 



(4.5) 



where t satisfies 



1 _ 1 - t , t{n - 2) 
P 



2 2n 

From this inequality and (|4.5j) . it follows that 



n(5-2) 

= c(i + ki|-V- 



Now, using this inequality along with 



and (|4.1j) . we obtain 



|V^| P < C 
< C 



+ Kl) f + Kl s (i + Kl 2 



2p + ,i(p-2) p 
L^O 4 | 2 I 1 



This completes the proof for n > 3. 

If n = 2, then from Sobolev's inequality, Holder's inequality, and (|4.4p . we have 



< 



< 



c 


(2 






C 


(2 


- ?)* 




c 



|V0|« dy 
\Vcb\ 2 dy) \n e 



(2-?) 



2 V 
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= 2±. 

9-2 

inequality, we obtain 



where q* — is the Sobolev conjugate of q. Then, again applying Holder's 



G 



(2-?) 
C 



(2-9)- 

Now using (|4~3)) , (g3), and (|4~Tj) . we obtain 

C 



(i + kil J 



|V<£| P < 



< 



(2-9)" 
C 



1 



(2-9)- 



qp + 2(p-q) 


p 




K\ 2 " H 


-Kl 5 - 





□ 



Lemma 4.2. for i/ie eigenfunction (j) e k , J < k < J+mj — 1, and any w 6 Hq(Qq), 
we have the following estimate: 



n(p-2) 



(4.6) 



where \\w\\i is from Lemma \jA\ with q(f, g) = / a^f^g^ dy, andC only depends 

Jn 

on \tto\, n, o~j, and is independent of e. 



Proof. First, recall that w is extended to be outside flo and 0| is extended to 
be in (B £ U B e ) n Vt° e . We have 



< 



affmnvew)? dy-o-i f {4>l) a {^w) a dy 



= \I + H\ 

+ \III + IV\ . 

First, since (f> £ k is an eigenfunction with eigenvalue er|, we have that III + IV = 0. 
Also, by Holder's inequality and Poincare's inequality, we have 



\i + n\<j \\m\mB sU B £ )(\\Vwi 



c 

_ '.m\\ L 2 (B e UB e ) 2 (B E UB E 

<c , l|v^|| z3(B , uge )Hli 



| W || 



I 2 (B t UB e ) , 
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where we have used Garding's inequality (|2.5[) on the last line for w. Thus, from 
Holder's inequality and Proposition 14. 2i 



\I + II\ < Ce^WVcPlh^Wwh 

Since <r° = Cj> the proof of the lemma is concluded. □ 

If we choose an interval / around such that af. € /, and let g) = 

a Hf f?9j dy and / = ?7 e 0|, we aim to satisfy the hypotheses for part 1 of 

Lemma 14.11 In order to do this, we need ||r7 e 0^ ||i 2 (o ) to be bounded away from 
0. To achieve this, we start with the following well-known proposition. 

Proposition 4.3. If A is an N x N matrix and v is a N x 1 vector such that 

N 

Av = and \Au\ < \Au\ for each I — 1, N, then v = 0. 

The next proposition shows that the functions {r]s4>k}k=T J ~ 1 are auTlos t or- 
thonormal. 

Proposition 4.4. For any e > and I, k £ N, (J < l,k < J+mj — 1), i/^fc = </>%, 
we have the following estimates: 

V 2 £ \M 2 dy^l-Ce^ (4.7) 



n 2 e (j) k ■ dy 

where C only depends on \£Iq\, n, and a®, and is independent of e 



d(p—2) 

<Ce^T~ ifk^l (4.8) 



Proof. We start by showing (|4.7|) . Since the eigenfunctions are normalized, we 
obtain for each k, 



V 2 M 2 dy= I (l-r? e )\<t> k \ 2 dy 



(1-^)|^| 2 dy 

T e UB e UB e 



, 0-3 



< \\VM\h(nJ T ^B £ UB £ \^T- 

< C k e p 

where, from (|4.2p . C k depends on <r°. Again, since a\ = Oj, we have (|4.7[) . 



1G 



Next, to show (|4.8p . we have 



r) e <t> k ■ 91 dy 



< 



V e 9k ■ 9i dy 

B E UB E 



rj e (j) k ■ (f>i dy 



n \(B E uB E 



Ve9k • 9i dy 



B E UB E 



Ok ■ 9l 



dy 



< I \9k ■ 9i\ dy + 

B E UB E J T E UB e UB. 



n \{B E uB E ) 

\4>k ■ <t>i\ dy, 



»k • 9i 



dy 



the second inequality following since the set of eigenfunctions form an orthogonal 
set in L 2 (fl £ ). So, next by Holder's inequality, we get 



r]l<l> k ■ 91 dy 



< 



\9k\ dy 

B E UB E 



dy 



B E UB E 



\9k\ 2 dy 



T E UB E UB E 



= 1 + 11. 

Now, from Poincare's inequality and (|4.2[) , we get 



\9i? dy 



T E UB E UB E 



I < 



dy 



B E UB E 



I-Br u b f 



dy 



B E UB E 



I -Bp U B F 



< \\V9k\\ L p(n E )£ 2f ||V0i|| i F (n , ) e ^ 

n(p-2) ^(p-2) 

< C k e 2 " Cie 2 p 

where C k again depends on a® and Ci depends on of. Thus, we have 



I < Ce—f— 

where C depends only on |f2o|, n, and a®. Similarly, 

d(p-2) 



II < Ce 



so that the proposition is proved. 



(4.9) 

(4.10) 
□ 



Note that with the aid of Lemma I4T2I and Proposition 14. 4[ if s is small enough, 
we have satisfied the hypotheses for part 1 of Lemma T4. II with q(f,g) = 
Iq a if f?9j dy an d / = T]e9%- Here, we relabel e% to be small enough to achieve 
this for any £ < £i, and note that E\ only depends on fixed parameters. To satisfy 
the hypotheses for part 2 of Lemma |4. 11 we need an orthonormal basis. The next 
proposition shows that for small e, we have a basis. 

Proposition 4.5. The set {'He4 ) k]k=j forms a linearly independent set for any 
N > J, for < e < Sq(N), where Eq(N) depends on N. 
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Proof. Assume Cjr/ £ (f>j + ... + CnVs^n ~ 0- Then, multiplying this equation by 
r/^f, we obtain 

N 
k=J 

So, if Aki = (?7e0|, r]e<t>f )L 2 (fi e )j we obtain by (14. 7p and 14.8(1 that 

d(p-2) 
d(p-2) 

> Ce » 

N 



if e < e(A), where e(A) depends on TV due to applying ([4.8(1 N — J times. Thus, 
we may use Proposition 14.31 to see that by setting C = (Cj, Cat)*, we have 
C = 0, so that the proposition is proved. □ 

Now we define J : L 2 (n e ) -> L 2 (ft ) to be given by Jo/ = ?7 £ /, and similarly, 
we define J e : L 2 (Qq) —> L 2 {Q, £ ) to be such that 



Jef(x) 



f(x), if x G fi 

0, if x e n e \fio- 



Let / = (a°j - Me^V 1 , aj+a J+™j j f or Af > to be chosen later. Also, let n 

be the projector onto the space spanned by the eigenfunctions corresponding to 
the eigenvalues, {crf^iLj, in /. We first consider e = e%. By Proposition 14.11 we 
may choose M = M(si) so that cr| is in / for J < fc < N , where N > J + mj — 1, 
and where N depends on e\. We next note that as e gets smaller, we may choose 
M = M{e) so that the set of eigenvalues in /, {<t|.}^j, will have index No in the 
range J + mj — 1 < Nq < N since our family {f2 £ } is nested. Our aim is to show 
that for e small, Nq = J + raj — 1. 

We apply Proposition ^. 51 to get the existence of £o(A) < £i so that {rje4'k}k=J 
is a linearly independent set for e < So(N) and for any N in the range J+mj — 1 < 
A < AT. Then, we choose M = M{s(N)) so that {r? e ^|}£j is also a basis for 
the range of JqTIJ £ . Thus, we may apply the Gram-Schmidt process to this basis. 
That is. define 



fj = V: 



eVj 



, _ , e (Ve4>%Jj) t (Ve(f>%, fk-l) t 

Jk — VePk II f i|2 J-l ■■■ iTf ill J k - 

\\fj\\ \\h-i\\ 



We have the following lemma: 
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Lemma 4.3. Let I be as defined above. For each k, J < k < J + mj — 1, we have 

d(p-2) 

Ce 4 ? 

||IL;c(/fc)||i < — , for e < e(N), and where M only depends on Oj, cr J _ 1 , 

and e(N). 

Proof. Following the previous arguments, when e = e(N), we find M = M(e(N)) 
so that {%</>|}^j is a basis for the range of JoTlJ £ , and then apply the Gram- 
Schmidt process to this basis. We note the dependence on a® and Oj_i is so 
that we only have 1 eigenvalue (with respect to fio) in /. So, defining q(f,g) — 
In a if ' f?9j dy, we may apply Lemma l4~2l Proposition 14. 4i and then Lemma |4~T1 
(part 1) to obtain 

d(p-2) 

Ce(N)^i^ 

» n '-<«»' £ -immr 

where C depends on \ flo\, n, crj}, and Oj +m ■ Then, from Proposition 14. 4[ Lemma 
14. 21 and properties of the norm, we get the result for e(N) and J < k < J+mj — 1. 
Then, for e < s(N), we may repeat this argument to get the result with e(N) 
replaced with e and M(e(N) replaced with M(e). But, since M(e(N)) < M(e), 
we obtain the desired result for e < e(N). □ 

We now let E = span{0|}^t™ J Also, let 11/ be the spectral projector 
corresponding to the eigenvalue <jj and rig be the spectral projector onto E. 

d(p-2) 

Ce 4 p 

Corollary 4.1. We have HII, - J U E J E \\ C {L^n„)} < , for e < e(N), 

where M only depends on Oj, 0j_i, and e{N). Consequently, for some e(J), 
No = J + m,j — 1 when e < e(J). 

Proof. Again, we first show for e — e(N). Normalize the ft-'s and observe that 

< d(--2) • Then apply Lemma 14.11 (part 2) to the normalized 

WfkW l~Ce(N)^r~ 

functions. Then for general e < e(N), we note that since Lemma T4.3I is true with 

d(p-2) 

Ce 4 p 

a uniform M, we obtain ||ILj — JoUe Je\\c{L 2 (n )} fs • We next note that 

NL 

if Nq > J + mj — 1 for all e < e(N), then we may find another projector LT^ 

d(p-2) 

Ce 4 p 

so that 1 1 11/ — Jo^-AJe\\c{L 2 (n )} < jj ■ But this would mean || JoH-E-h — 

d(p-2) 

J IUJe||.c{L 2 (r2o)} ^ ■ Therefore, for some e(J), N — J + mj - 1 when 

e<e{J). ' □ 



Proof of Theorem \2.1\ We first prove for J = 1. By Corollary [4TTJ for e < e(l), we 

d(p-2) 

obtain mi = Nq. This implies that |af — cr^| < Ce 4 p only for k, 1 < k < mi, 
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and hence, the result for J = 1. The result for J = 1 implies that not only may we 
choose M so that all eigenvalues {&k}™=^+i are in the interval corresponding to 
the next highest eigenvalue a^ ni+1 , but also that a® is not in this interval. Thus, we 
apply the same reasoning here to get the result for cr^ i+1 . Then, by an induction 
argument, we get the result for each J e N, satisfying <Tj > <Tj_ 1 . We note here 
that since C depends on e(J), it depends on the multiplicity J. □ 

We note that this paper introduces the use of L p -estimates obtained by the 
reverse Holder technique to the study of spectral problems for elliptic operators. 
Thus, this technique may be useful in studying spectral problems in situations 
where we do not know if higher regularity of solutions is true. We close by listing 
some open problems. 

• If we have some additional regularity on the domain, can we use the methods 
from this work to get convergence of Neumann eigenvalues for general elliptic 
systems? 

• For elliptic systems on a symmetric dumbbell region with a straight tube, 
can we achieve upper and lower bounds on the splitting between the smallest 
eigenvalues? 

• Can we investigate this problem further to see if a better rate of convergence 
exists? 
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